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Abstract 

A brief introduction is given in the generic microscopic framework of superconductivity. The 
consequences for the temperature dependence of the kinetic energy, and the correlation energy 
are discussed for two cases: The BCS scenario and the non-Fermi hquid scenario. A quantitative 
comparison is made between the BCS-prediction for d-wave pairing in a band with nearest neighbor 
and next-nearest neighbor hoppping and the experimental specific heat and the optical intraband 
spectral weight along the plane. We show that the BCS-prediction produces the wrong sign for the 
kink at Tc of the intraband spectral weight, even though the model calculation agrees well with 
the specific heat. 
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FIG. 1: Experimental internal energy and the specific heat obtained by Loram et azQ]. The 
extrapolations are obtained following the procedure of Ref. 

I. MODEL INDEPENDENT PROPERTIES OF THE SUPERCONDUCTING 
STATE 



A. Internal energy 

When we cool down a superconductor below the critical temperature, the material enters 
a qualitatively different state of matter, manifested by quantum coherence over macroscopic 
distances. Because the critical temperature represents a special point in the evolution of 
the internal energy versus temperature, the internal energy departs from the temperature 
dependence seen in the normal state when superconductivity occurs. Because for T < Tc 
the superconducting state is the stable equilibrium state, the internal energy in equilibrium 
at T = is an absolute minimum. Hence cooling down from above the phase transition 
one would expect a drop in the internal energy when at the critical temperature. This 
drop of internal energy stabilizes the superconducting phase among all alternative states of 
matter. An experimental example of this well known behavior is displayed in Fig. Q where 
the internal energy was calculated from the electronic specific heat data|]]] according to the 
relation EintiT) = c{T')dT'. The broadened appearance of the phase transition suggests 
that the superconducting correlations disappear rather gradually when the temperature is 
increased above the phase transition, an effect which remains noticeable in the specific heat 
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graph up to 125 K. 

Understanding the mechanism of superconductivity means to understand what stabihzes the 
internal energy of the superconducting state. In BCS theory superconductivity arises as a 
result of a net attractive interaction between the quasi-particles of the normal state. Note, 
that implicitly this approach is firmly rooted in the paradigm of a Fermi-liquid type normal 
state. On the other hand, the school based on Anderson's original work|2] asserts, that a 
strong on-site repulsive interaction can also give rise to high Tc superconductivity. The latter 
models typically require that the material is not a Fermi liquid when superconductivity is 
muted, either by raising the temperature or by other means. 

B. Pair-correlations 

Without loss of generality, i. e. independent of the details of the mechanism which leads 
to superconductivity, it is possible to provide a microscopic definition of the superconducting 
state. For this purpose let us consider the correlation function defined as 

G{r,R,,R,) = (V'}(i?i+r)^|(i?0^i(WT(^2 + r)) (1) 

where ipKRj) is a single electron creation operator. This function defines two types of 
correlation: (i) the electron pair-correlation as a function of the relative coordinate r, and 
(ii)v the correlation between two pairs located at center of mass coordinates Ri and i?2- 
In the normal state there is no correlation of the phase of G{r, Ri, R2) over long distances 
\Ri — -R2I due to the finite mean free path of the electrons. As a result the integral over the 
center of mass coordinates of the correlation function 

= ^/ "^'^1 / d'RMr,Ri,R2) (2) 

averages to zero in the normal state. In contrast the superconducting state is characterized 
by long range phase coherence of the center of mass coordinates, implying (among other 
things) that the correlation function averaged over all center of mass coordinates, g{r), is a 
finite number. 
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II. INTERNAL ENERGY AND ITS DECOMPOSITION USING THE BCS 
MODEL 



A. BCS: Correlation function, and correlation energy for d-wave pairing 

In the weak coupling scenario of BCS theory the electrons have an effective attractive 
interaction, as a result of which they tend to form pairs. For the purpose of the present 
discussion we will assume that the interaction is of the form 

W = V{r - r')n{r)h{r') (3) 

where n{r) is the electron density operator. The interaction energy in the superconducting 
state becomes lower than in the normal state, due to the fact that the effective attractive 
interaction favors a state with enhanced pair-correlations. The value of the interaction 
energy of the superconducting state, relative to the normal state is 

(H^)^ - (W)^ = I dhg{r)V{r) = Y,9kVk (4) 

k 

where Qk and Vk are the Fourier transforms of g{r) and V{r) respectively. Using the Bogoli- 
ubov transformation the correlation function can be expressed j J] in terms of the gap-function 
Ak and the single particle energies Ek = {{ek — f^Y + A|}^/^. 

This corresponds to the conversion of a pair (g, —q) to a pair with quantum numbers 
[q + k, —q — k). In the expression for the correlation energy, Eq. |3J the transferred mo- 
mentum k is carried by interaction kernel Vk- 

Starting from a model expression for the single electron energy-momentum dispersion e^, 
and the gap-function A^, it is a straightforward numerical exercise to calculate the summa- 
tions in Eq. Adopting the nearest neighbor tight-binding model with a d-wave gap, and 
adopting the ratio A(7r, 0)/iy = 0.2, where W is the bandwidth, the correlation function 
Qk can be easily calculated, and the result is shown in Fig. |21 We see from this graph that 
a negative value of {H^) ^ — (-f^*)„ requires either (i) < for k in the neighborhood of 
the origin, or (ii) > for k in the vicinity of (7r,7r). The corresponding representation 
in real space, g{r), shown in Fig. |2l illustrates that the dominant correlation of the d-wave 
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FIG. 2: The k-space (top panel) and coordinate space (bottom panel) representation of the su- 
perconductivity induced change of pair-correlation function for d-wave symmetry (bottom panel). 
Parameters: A/VF = 0.2, wd/VF = 0.2. Doping level: x = 0.25 

superconducting state is of pairs where the two electrons occupy a nearest neighboring site, 
while the on-site amplitude is zero. 

Combining the information of Fig. |21 with Eq. |2l it is clear that the strongest saving of 
correlation energy is expected if the electrons interact with an interaction of the form 
V{ri, = Vq J2a "^(^1 — ^2 + a) where the vector a runs over nearest neighbor sites, and Vq 
is a negative (meaning attractive) interaction between electrons on nearest neighbor sites. 

B. BCS: The gap-equation, specific heat and internal energy 

To illustrate the predictions of BCS theory for the temperature dependence of the corre- 
lation energy and the kinetic energy, we start by solving the gap equation 



which must solved together with the constraint that the average number of particles is tem- 
perature independent. This requires that the chemical potential must be adjusted to keep 
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FIG. 3: BCS prediction of the d-wave gap-function and the chemical potential. 



the thermal average of Y^/^ {hk) = N aX a, constant value. In the numerical examples of this 
paper we have adopted N = 0.85 corresponding to a hole doping of 0.15. This solution of the 
gap equation is shown in Fig. El Notice, that the temperature dependence of the chemical 
potential follows closely the experimental observations reported for Y123 in Ref. [3| 
To check that the band-parameters used here are reasonable, we display in Fig. |3] the corre- 
sponding prediction for the specific heat and the internal energy, using the same parameters 
as for Fig. El If we compare this to Fig. HI we see that the band-parameters adopted 
here quantitatively reproduce the observed specific heat. Hence the present set of band- 
parameters {t,t') and the coupling parameter Vq represent the best phenomenological choice 
for quantitative testing of the BCS-model. 

C. BCS: Temperature dependence of the correlation energy and the kinetic energy 

The BCS prediction for the temperature dependence of the average interaction energy 
follows from Eq. We then use the BCS variational wave-function for the statistical average 
of Eq. resulting in 
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Input parameters for BCS model: 
^kq = ^0 (cosk^ - cosk )(cosq^ - cosq ) 

= tcosk^ + tcosk + 2t'cosk^cosk 
V„ = 88 (r^ieV) 
t = -297.6 (meV) 
f= 81.8 (meV) 
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FIG. 4: BCS prediction of the internal energy and the specific heat. 



Simultaneously there is an increase of the 'kinetic' energy 



tanh 



Ek 



2kBT 



(8) 



In Fig. this is displayed, using the same parameters as for Fig. |21 



III. RELATIONSHIP BETWEEN INTRA-BAND SPECTRAL WEIGHT AND KI 
NETIC ENERGY 

A measure of the kinetic energy is provided by the following relationj^ Q| 



/ dujRea{uj)duj = vr— — 22i^k,a) 



(9) 



where the high frequency limit indicates that the integral should include only the intra- 
valence band transitions, and the condensate peak at = if the material is a supercon- 
ductor. The integral over negative and positive frequencies (note that (7(00) = o'*{—uj)) avoids 
ambiguity about the way the spectral weight in the condensate peak should be counted. If 
the band structure is described by a nearest neighbor tight-binding model, Eq. El leads to 
the simple relation 

PL = j^^ Rea{u)duj = {-Hk^n) (10) 
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FIG. 5: BCS prediction of the kinetic energy and the correlation energy. 

Hence in the nearest neighbor tight-binding Hmit the partial f-snm provides the kinetic energy 
contribution, which depends both on the number of particles and the hopping parameter 
tjlll. Ql- However, if the band-structure has both nearest neighbor hopping and 

next nearest neighbor hopping, Eq. El is not an exact relation, and instead Eq. IHl should be 
compared directly to the experiments. In Fig. IHlwe compare the spectral weight, calculated 
directly using Eq. El to the result of Eq. [TUl using the same parameters as for Fig. 01 
Note that the kinetic energy has to be divided by a factor two, as we are interested in the 
projection along one of the two axes in the ab-plane, which can be compared directly to the 
experimental value of pl- From Fig. (HI we can conclude, that the effect of including t' in the 
calculation is rather small, and it is still OK to identify pl with the kinetic energy apart 
from a minus sign. 



IV. EXPERIMENTAL DETERMINATION OF THE INTRA-BAND SPECTRAL 
WEIGHT 



In two recent experimental papers measurements of pi in Bi2212 have been reported 
The values of the kinetic energy change in the superconducting state were in 
quantitative agreement with each other, and both papers arrived at the same conclusion: 
Contrary to the BCS prediction, the kinetic energy of the superconducting state is lower 
than in the normal state (taking into account a correction for the temperature trends of the 
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FIG. 6: BCS prediction of the spectral weight function. 

normal state). In Fig. [3 the data of Ref. j3] have been reproduced. Comparing this with 
the BCS prediction clearly demonstrates the large qualitative discrepancy between theory 
and experiment. Clearly the type of mechanism assumed in BCS theory is not at work here. 



V. IMPLICATIONS OF THE EXPERIMENTAL DATA 

The trend seen in the experimental data has been predicted by Hirsch in 1992^, 18. igj]. 
The model assumption made by Hirsch was, that the hopping probability of a single hole 
between two sites becomes larger if one of the two sites is already occupied by a hole. 
Although this model provides good qualitative agreement with the optical experiments, it 
has one serious deficiency: It also predicts s-wave symmetry for the order parameter, in 
sharp contrast to a large body of experimental data which show that the superconducting 
gap in the cuprates has d-wave symmetry. 

In a recent set of calculations based on the Hubbard model, Jarrell et aljs^ obtained a 
similar effect as seen in our experiments, both for underdoped and optimally doped samples. 
Crudely speaking the mechanism is believed due to the frustrated motion of single carriers 
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FIG. 7: Experimental values of the ab-plane spectral weight function, taken from Ref. 

in a background with short-range (RVB-type) spin-correlations, which is released once pairs 
are formed. 

VI. CONCLUSIONS 

We have made a quantitative comparison between the BCS-prediction for d-wave pairing 
in a band with nearest neighbor and next-nearest neighbor hoppping and various experi- 
ments, in particular specific heat and measurements of the optical ab-plane sumrule. We 
have shown that the BCS-prediction produces the wrong sign for the kink at Tc of the ab- 
plane intraband spectral weight, while the model calculation is in good agreement with the 
experimental specific heat data. 
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